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Abstract 

An approach is proposed to calculate the direct reaction (DR) and fusion 
probabilities for heavy ion collisions at near-Coulomb-barrier energies as func- 
tions of the distance of closest approach D within the framework of the optical 
model that introduces two types of imaginary potentials, DR and fusion. The 
probabilities are calculated by using partial DR and fusion cross sections, 
together with the classical relations associated with the Coulomb trajectory. 
Such an approach makes it possible to analyze the data for angular distribu- 
tions of the inclusive DR cross section, facilitating the determination of the 
radius parameters of the imaginary DR potential in a less ambiguous manner. 
Simultaneous x 2 -analyses are performed of relevant data for the 16 O+ 208 Pb 
system near the Coulomb-barrier energy. 
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I. INTRODUCTION 



Collisions between heavy ions at near- Coulomb-barrier energies are very much governed 
by the Coulomb potential involved, and thus the general features of the elastic scattering 
and direct reaction (DR) data can be understood based on the idea that the colliding ions 
primarily move along a classical Coulomb trajectory Jl|. These features are seen most 
dramatically in plots of the ratios of the elastic differential cross section (daE/dQ) or the 
inclusive (sum of all different) DR one (dao/dfl) to the Rutherford differential cross section 
(da c /dfl), i.e., 

« S Sf^ <= « = EorD), (1.1) 

as a function of the distance of closest approach D (or the reduced distance d) that is 
related to the scattering angle 9 through 

D = d{A^ + Ar) = \DJl + — )—) with D = ^. (1.2) 

2 \ sin(0/2)/ E cm 

Here D is the distance of closest approach in a head-on collision (s wave). Further, (Ax, Z{) 

and (A 2 , Z 2 ) are the mass and charge of the projectile and target ions, respectively, and E cm 

(Ei a b) is the incident energy in the center of mass (laboratory) system. Pe and Pd thus 

defined may be called the elastic and DR probabilities, respectively. The impact parameter 

b and orbital angular momentum I, specifying the trajectory, are related to 9 and D by 

b = I = ^ot°- = sjD(D-D ), (1.3) 

where k is the wave number. 

As an illustration, we show in Fig. 1 such plots of Pe and Pd for the 16 O+ 208 Pb sys- 
tem HQ at five different incident energies of .5^=80, 84, 90, 96, and 102 MeV for P E and 
at a single energy of Ei ab =90 MeV for P D , where the data are available. As seen, the values 
of Pe at different energies line up to form a very narrow band and take a value very close 
to unity for, say, d > 1.65 fm (= dj, interaction distance). When d becomes smaller than 
dj, Pe falls off very rapidly, approximately exponentially. 



The observed behaviour of Pe may easily be explained based on the physical picture 
that the projectile ion moves primarily along a Coulomb trajectory. For the case d > dj, the 
trajectory is far away from the target and the projectile is scattered at the Coulomb scat- 
tering angle 9 without being influenced at all by the nuclear force. The resultant scattering 
cross section is thus equal to the Rutherford cross section, and Pe becomes unity. When d 
becomes smaller than di, however, the incident ion gets under the influence of the strong 
nuclear interaction, and absorption takes place, reducing the P^-value below unity. 

In accordance with the observed behaviour of Pe, Pd starts to have a significant value 
at d ~ dj. It reaches its maximum value (Pd ~ 0.24) at d ~ 1.58 fm, where Pe becomes 
approximately 0.7. In the region of d = 1.58 ~ 1.65 fm, the sum Pe + Pd thus stays close 
to unity. This indicates that in that region, the main cause of absorption in the elastic 
channel is DR processes. When d becomes still smaller, the sum Pe + Pd falls off rapidly 
from unity, showing that absorption due to more complicated processes eventually leading 
to fusion takes place. It is remarkable that the sum Pe + Pd becomes extremely small, say, 
10~ 3 and thus essentially zero, at around d = d c = 1.30 fm, which is the distance for the s- 
wave Coulomb-barrier top. This means that the incident flux is almost completely absorbed 
when it reaches that distance. It is worth noting that the same picture holds irrespective of 
the incident energy, so long as it is not far away from the Coulomb-barrier energy. 

Theoretically, we have a very well established optical model for evaluating Pe- This is 
not the case for Pd- There are a variety of theoretical methods that has been proposed for 
calculating contributions from inelastic scattering and transfer reaction processes to Pd by 
means of either semi-classical or classical approximations Jl[]. It is, however, still a formidable 
task to carry out calculations including all possible processes to obtain a theoretical value of 
Pd- The aim of the present work is to propose a simple approach to calculate Pd within the 
framework of an optical model that introduces two types of the imaginary potentials; one 
for DR, the other fusion 0-0. We propose to evaluate Pd from the partial absorptive DR 
cross sections generated from the optical model calculation. The underlying assumption is 
that even after the reaction (removal from the elastic channel) the projectile ion still moves 
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along the Coulomb trajectory, being eventually emitted at the Coulomb scattering angle. 
Under the assumption, we may use the classical relation Eq. (1.3) to convert the partial 
wave cross section to Pp. 

The same prescription may also apply to the partial fusion cross sections. This enables us 
to evaluate the hypothetical fusion probability P F just as we calculate Pp. The conventional 
wisdom assumes [§] that 

P E + P D + P F ^ 1, (1.4) 

which expresses a simple physical idea that what is absorbed in the elastic channel at r = D 
goes into either DR or fusion channels. Since there is no such measured Pp-value available, 
it has been impossible to test the above relation experimentally. However, it is possible to 
examine the relation by using the theoretical Pp. In the present work, theoretical expressions 
for Pp and Pp are derived in Sec. II, where we also perform numerical calculations of Pp 
and Pp along with Pp and examine the validity of the relation Eq. (1.4). The feasibility 
of evaluating Pp enables us to analyze the angular distribution of the inclusive DR cross 
section. This facilitates the determination of the DR part of the optical potential with less 
ambiguity. We shall demonstrate this also in Sec. II. In Sec. Ill, we repeat the simultaneous 
X 2 -analyses that we made several years ago for the data on the 16 O+ 208 Pb system shown 
in Fig. 1. The reanalyses are needed since the fusion data have been revised || after Ref. @ 
was published. Sec. IV will then be devoted to our conclusions. 

II. OPTICAL MODEL CALCULATIONS OF P D AND P F 

A. Derivation of Theoretical Expressions for Pp and Pp 

In this subsection, we try to derive theoretical expressions for Pp and P F within the 
framework of the optical model. We follow the approach proposed some time ago to calculate 
the total DR and fusion cross sections within the optical model by using imaginary, surface 
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type DR and volume type fusion, potentials, Wo(r) and Wp(r), respectively f5|-[7|. The 
total DR and fusion cross sections are then calculated as PJTOfl 

(Ti = ^ <X (+ V,(Olx (+) > (t = DoiF), (2.1) 

where x is the usual distorted wave function that satisfies the Schrodinger equation with 
the full optical model potential U(r). <jp> and ap are thus calculated within the same 
framework as the differential elastic scattering cross section, dap/dfl, is calculated. Such a 
unified description enables us to describe all different types of reactions on the same footing. 

The basic ingredients for obtaining theoretical expressions for Pp> and Pp are the partial 
wave cross sections, a^i {% = D or F), which are obtained by simply expanding the cross 
sections, Eq. (2.1), into the partial wave components; <7j = Y^e a r,e- °V can explicitly be 
given as || 

a iit = ^(2£ + 1)T M {i = D or F), (2.2) 

where 

T i; ,= — / Mr)\ 2 Wi{r)dr. (2.3) 
hv Jo 

In the above equation, xe( r ) i s the partial distorted wave function and v is the relative 
velocity. 

Eq. (2.2) with Eq. (2.3) is still a quantum mechanical expression, where i takes only 
integer values. In what follows, we introduce a few semi-classical approximations customarily 
used |I|P, pJ]jr2]| . The first is to treat t as a continuous variable and to assume that 

^ = a, t . (2.4) 

We then use the classical relation Eq. (1.3) that relates i to 9. It is then straightforward to 
get 

da,i(i) 1 d£da t (£) kD 1 



dfl 27ism8d8 dt 16tt cos(fl/2)sin 3 (fl/2) 



CTi;£. (2.5) 



Inserting further Eq.(2.2) into Eq.(2.5), and dividing the resultant expression by the Ruther- 
ford cross section a c , one finally obtains 

21 + 1 

Pi = -rjr- tan(0/2)% « T i;£ (2.6) 

In obtaining the last expression, use is made of the approximation that 21 + 1 f» 2£. 
In order that Pi can be a probability, it should satisfy 

P 4 < 1. (2.7) 

This requirement is indeed satisfied; in fact we have 

P D + P F = T D;i + Tp.£ = T e = 1 - |^| 2 , (2.8) 

where Tg is the transmission factor and Sg is the partial wave 5-matrix. Since both Pp and 
Pp are positive quantities, it is clear from the above relation that Pp and Pp should be less 
than unity. (Note that there is no reason that Pp should be smaller than unity. Quantum 
effects such as interference and diffraction may cause the value to be greater than unity.) 
Now, for very small I- values, we expect that Pp, — > 0, hence 

P F ml- \S e \ 2 -> 1 for small i. (2.9) 

The last relation follows from the fact that for such a strong absorptive case as in heavy 
ion collisions, St becomes essentially zero for small £. Since Pp + Pp — > for small £, 
Pp + Pp, + Pp — > 1 as expected earlier in Eq. (1.4). In the next subsection, we further study 
this point numerically. 

In passing, we remark that the procedure we have proposed can also be used to reduce 
the quantum mechanical Rutherford cross section to the classical one. As is well known, 



this reduction has been given by using a set of semi-classical approximations [pTT| , [12|] . The 
quantum mechanical Rutherford cross section has the well known form, 
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-L^(2£ + l)( e 2 ^-l)P^) 



(2.10) 



where T)n is the Coulomb phase shift and Pi is the Legendre function. One of the approxima- 
tions introduced in the reduction process is to ignore the term — X^(2^ + l)Pe(9). This 
term gives rise to a divergent contribution at extremely forward angles and we ignore it as 
is usually done (see Refs. fTtUT^! ). We then integrate Eq. (2.10) over angles to obtain the 
total elastic cross section expressed as a sum of the partial cross sections, which is in turn 
converted to an integral over I. The resultant total elastic scattering cross section takes a 
very simple form, namely 

a c = £ a c , = J ^di with ^ = a c , = £(2£ + 1) « £(2*). (2.11) 

By inserting the last expression in Eq. (2.11) into Eq. (2.5), we obtain the Rutherford 
differential cross section. 

It appears that the procedure used for reducing the quantum mechanical Rutherford 
cross section to the classical one involves a contradictory element; we first integrate the 
differential cross section over angle, but then recover it from the partial wave cross sections 
we obtained as a result of the angle integration. However, the procedure can be justified: 
The quantum mechanical cross section Eq. (2.10) is given as a coherent sum over t. As 
has been demonstrated in a number of semi-classical treatments of Eq. (2.10) |TT| , |T2"fl , the 
dominant contribution to the differential cross section for a given scattering angle 9 comes 
from the partial waves around i = £g, where is related to 9 by Eq. (1.3). The contribution 
becomes 5-function-like in the classical limit of h —>■ 0. In the present procedure, we carry 
out an integration over 9 first, but from what has been discussed above, it is seen that the 
contribution from the angle 9 is stored into the partial wave cross section of I ~ Iq. It is 
thus justifiable to recover the differential cross section at an angle 9 from the partial wave 
cross section for t = £g. This procedure is most justified when the scattering is closest to 
classical. 
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B. Numerical Examples 



Calculations of Pp, Pd, an d Pp are performed for the 16 O+ 208 Pb system with incident 
energy £^=90 MeV, using the optical model potential as fixed in our previous study J7|. We 
present the results in Fig. 2, where the solid, (thick and thin) dotted, and (thick and thin) 
dashed curves are the calculated values of Pe, Pd an d Pf, respectively The experimental 
data for Pe and Pd are also plotted by the solid and open circles, respectively. As seen, 
the calculated Pe reproduces the experimental Pe very well. This is not the case, however, 
for Pjj\ the calculated Pd is shifted to the smaller <i-region by about 0.05 fm as compared 
with the experimental data, particularly in the region of d > 1.6 fm. Thus, the comparison 
of the calculated Pd with the data provides an additional test of the parameters used in 
Ref. @. In fact, this shift means the radius parameter td = 1-50 fm used in Ref. [|] is too 
small to describe the data. We thus repeated the calculation with a larger radius parameter 
of r D = 1-55 fm. The calculated Pd thus obtained is plotted by the thin dotted curve 
shown in Fig. 2. The fit to the data is improved. It is worth noting that the recalculated 
Pd shifts toward the larger d-region by 0.05 fm, the same amount as the increase of the 
radius parameter td- This shows that the experimental Pd provides a very sensitive test 
of the revalue. Based on the result obtained above, we use the value rp>=1.55 fm in the 
X 2 -analysis discussed in the next section. 

In the calculation of Pp shown in Fig. 2, use is made of rp=1.40 fm. Thus the Pp 
curve lies in much smaller <i-region than the Pd- We also observe that the slope of Pp is 
much steeper than that of Pd- This reflects the fact that the diffuseness parameter used 
for Wp{r]E) (ap=0.25 fm) is smaller than that of Wn(r;E) (ap = 0.45 fm). To show the 
effects of the dp- value, we present by the thin dashed curve another Pp calculated with 
ap=0A5 fm. It is seen that the slope of Pp at large distances is almost the same as that of 
Pd with cid = 0.45 fm. 

Let us now turn to the sum Pe+Pd+Pf shown by the dash-dotted curve. As expected, it 
stays very close to unity, confirming that the relation Eq. (1.4) is fairly well satisfied, within 
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the accuracy of, say, 20 %. The sum shows some oscillations around unity, which may be 
ascribed to quantum interference effects. The oscillation is also visible in the experimental 
Pe values. Accumulation of more accurate data may enable us to test this explanation in a 
more detailed manner. 

Presented in Fig. 3 are the P^-values calculated for the incident energies considered in 
Fig. 1. Use is made of the optical potential determined from the x 2_ana ly s i s discussed in 
the next section. Since we use such a potential as determined from the x 2 -^-, the calculated 
Pe fit the data given in Fig. 1 very well and thus they forms a band very much similar to 
that seen in Fig. 1. 



III. x 2 - ANALYSES 

We have repeated simultaneous x 2 - ana ly ses as i n Ref. J7J for the elastic scattering, DR, 
and fusion data for the 16 O+ 208 Pb system. This is motivated for two reasons: The first 
is that the fusion data have been revised ||, after Ref. |7J was published. The second is 
that we are now able to test the calculations against the data for Pjj. As described in the 
previous section, the value of the radius parameter r\D=1.50 fm used in Ref. [0] is too small 
to explain the data. A better revalue is Tb =1.55 fm. Other parameters must be fixed 
with this more appropriate value of rp. As in Ref. [[/]], we utilize a dispersive type optical 
potential [|13[ 



U = U c {r) - [V {r) + V{r; E) + iW{r; E)}, (3.1) 

where Uc{r) is the Coulomb potential and Vo(r) is the energy independent Hartree-Fock 
part of the potential, while V(r; E) + iW(r; E) is the polarization part of the potential f!4 



that originates from couplings to reaction channels. They are assumed to have volume-type 
fusion and surface-derivative-type DR parts. Explicitly, Vo(r), V(r;E) and W(r;E) are 
given, respectively, by 

V (r) = V f(X ), (3.2) 
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V(r; E) = V F (r; E) + V D (r; E) = V F (E)f(X F ) + AV D (E)a D ^§^-, (3.3) 

W(r; E) = W F (r; E) + W D (r; E) = W F (E)f(X F ) + 4W D (E)a D ^^-, (3.4) 

where f(Xi) = [1 + exp(Xj)] -1 with Xi = (r — Rij/di (i — 0, D or F) is the usual 
Woods-Saxon function. Use is made of the values used in Ref. for the parameters of the 
bare potential Vo(r), and the geometrical parameters of V(r; E) and W(r;E) (except 
as discussed above); Vo=60.4 MeV, r =1.176 fm, a =0.658 fm, rp=1.40 fm, ap=0.25 fm, 
ro=1.55 fm, and aD=0.45 fm. Once the geometrical parameters are fixed, the dispersion 
relation is reduced to a relation for the strength parameters Vi(E) and Wi(E) (i = D and 
F). The relation now reads |U 



Vi{E) = Vi (E,) + r de- Zf2 57. (3-5) 

v ; v ; 7T Jo {E> - E S )(E' - E)' v ; 

where P stands for the principal value and Vi(E s ) is the potential value at a reference energy 
point E = E s . 

As was done in Ref. [0], we approximate the ^-dependence of Wi(E) just above the 
threshold energy E 0]i (defined as Wi(E 0]i ) = 0) by a linear function of E. (See the forthcom- 
ing Eqs. (3.7) and (3.8).) We then identify this threshold energy as that determined from 
the linear representation of the quantity S{(E) introduced by Stelson et al. [jHJ as 



Si = ^/Ea t (E) oc (E - E . t ) (z = D or F). (3.6) 

The threshold energies E .i thus defined are essentially the threshold energies of the DR 
(i = D) and fusion (i = F) cross sections, and it is plausible to identify the two threshold 
energies to be the same. The authors of Ref. [|l6j considered the quantity Si only for the 
i = F case, but we extend it to DR. Originally, two threshold phenomena in the imaginary 
part of the optical potential and the fusion cross section data were found independently, but 
it was noticed later [17]] that the two are very close to one another. Once we have separated 



the imaginary potential into the DR and fusion parts, it is physically plausible to require 
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that the two thresholds should be the same. In Fig. 4, we present the Si-values for i = D 
and F. There we find that E -p =73.0 MeV and E Q .p=76.0 MeV, which will be used later 
as the threshold energies of Wp(E) and Wf(E), respectively. 

In an attempt to determine the polarization potential, simultaneous x 2_ana ly ses were 
performed, treating all four strength parameters, Vp, Vp, Wp, and Wp as the adjustable 
parameters. We took into account all the data |3],^| available for incident energies between 
Ei a b=S0 MeV and 104 MeV. We included the total DR and fusion cross sections in the 
analyses. 

The values of the parameters thus extracted are presented in Fig. 5 for Vp and Wp, and 
in Fig. 6 for Vp and Wp. Let us consider first the results for Vp> and Wp,. A considerable 
fluctuation is seen in the values of Vp, but Wp, changes smoothly as a function of E. The fact 
that Wp could be fixed as a smooth function of E indicates that these values are reliable. 
There is a reason that Wp can be determined rather unambiguously and becomes a smooth 
function of E. It is because Wp is the dominant absorptive term in the peripheral region. 
Therefore, the elastic scattering cross section is quite sensitive to the value of Wp. This is 
not the case for Vp] at the strong absorption radius, where the elastic scattering cross section 
is sensitive to the real potential, Vp is generally much smaller than the bare potential Vo(r), 
resulting in difficulty in determining Vp unambiguously. The fluctuation seen in Fig. 5(a) 
may be understood to arise from this difficulty. 

The IVo-values determined from the x 2 -analyses can be well represented by the following 
function of E (in units of MeV): 

for E <73.0 

W D (E) = \ 0.015(^-73.0) for 73.0< E <92.5 (3.7) 
0.2925 for E > 92.5 

where _Eo ; d=73.0 MeV is used as extracted in Fig. 4. The solid line shown in Fig. 5(b) is 
Wp given by Eq. (3.7). The line fits the empirical values quite well. 

Since a reliable value of Wp is now available, one can calculate Vp by using the dispersion 
relation Eq. (3.5). In doing this, we need to know one more parameter, i.e., the value of 
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Vd at E — E s . We may fix this Vd(E s ) by fitting the average of the resultant Vp to that 
of the empirically determined Vp. The solid curve shown in Fig. 5(a) shows the Vd- values 
thus calculated. The V D {E s )-\ahxe used is V D (E a )=QA MeV at £ s =92.5 MeV. 

As seen in Fig. 6, Vf and Wf are both determined as fairly smooth functions of E. The 
Wp-vdlues may be represented (in units of MeV) as 

for E <76.0 

W F {E) = { 0.32(£ - 76.0) for 76.0< E <86.0 (3.8) 
3.2 for E > 86.0 . 

Again we took the threshold energy of E . F =76.0 MeV determined from Sp- The solid 
line shown in Fig. 6(b) represents Wf in Eq. (3.8). We then calculated the V F (E) by 
using the dispersion relation Eq. (3.5) with Wp given by Eq. (3.8). The reference potential 
V F {E S ) involved was chosen as Vf(E s )— 3.50 MeV at E s =8Q.O MeV. As shown by the solid 
curve in Fig. 6(a), the predicted Vp- values again agree reasonably well with the empirically 
determined values. 

We take as our final potential parameters Wp, and Wf given, respectively, by Eqs. (3.7) 
and (3.8), and also Vp and Vf generated from the dispersion relation Eq. (3.5). The potential 
with such parameters then fully satisfies the dispersion relation. Using such an optical 
potential, we calculated the final theoretical Pp, <jd, and <jf and presented them in Figs. 7 
and 8 in comparison with the experimental data. As seen, all experimental Pp, up and dp 
are well reproduced by the calculations. 

We now wish to give some remarks on the polarization potential we obtained. First, 
there is a remarkable difference in the the energy dependences between the DR and fusion 
potentials. A very rapid change is seen only in the fusion part of the potential. The slope 
of Wf(E) given by Eq. (3.8) in the threshold region is 0.32, while that in Wd(E) given by 
Eq. (3.7) is only 0.015. As a result, we see a significant energy variation of about 2 MeV 
in Vf(E) in the interval of ~ 10 MeV, but the change in Vp is only 0.1 MeV in the energy 
range of ~ 20 MeV. We may thus conclude that the threshold anomaly exists in the fusion 
part of the potential, but not in the DR part. 
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IV. CONCLUDING REMARKS 



We have presented a simple method to calculate the DR and fusion probabilities within 
the optical model by introducing two types of imaginary potentials, DR and fusion. These 
probabilities are calculated by using the partial DR and fusion cross sections generated from 
the corresponding imaginary potentials with the help of the classical relation between the 
orbital angular momentum t and the scattering angle 9. The probabilities thus calculated 
were shown to satisfy the condition that the value should be equal to or less than unity. 

Based on the expressions derived, numerical calculations of these probability were per- 
formed. We found that the sum of the DR, fusion and elastic probabilities stays close to 
unity. We also analyzed the angular distribution data of the inclusive DR cross section, 
demonstrating that the data provide useful information for determining the radius param- 
eters of the DR potential. It was observed that a very rapid energy variation (threshold 
anomaly) was in the fusion part of the potential, but it is hardly seen in the DR part, 
particularly in the real part of the potential. 

A simultaneous x 2 -a n alysis of elastic scattering, DR and fusion cross sections for the 
16 O+ 208 Pb system at near-barrier energies were performed for determining the polarization 
part of the optical potential that satisfies the dispersion relation over all space. The potential 
thus determined is found to reproduce the data well. 

The authors wish to express their sincere thanks to Prof. W. R. Coker for his kind 
reading of the manuscript and comments. One of the authors (BTK) acknowledges the 
support by Korea Research Foundation (KRF-2000-DP0085). 
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FIGURE CAPTIONS 



Fig. 1. The experimental elastic probabilities, P E , as a function of the reduced distance 
D for the 16 O+ 208 Pb system at E lab = 80, 84, 90, 96, and 102 MeV. The measured DR 
probabilities, Pp, are also plotted at E [ab = 90 MeV. The data are taken from Refs. 3 and 
4. 

Fig. 2. The optical model predictions of elastic, DR, and fusion probabilities as a func- 
tion of the reduced distance D for the 16 O+ 208 Pb system at E iah = 90 MeV are shown in 
comparison with the experimental ones. The theoretical total probabilities are also shown. 
The thick dotted curve denotes Pd calculated with the dispersive optical potential deter- 
mined in Ref. 7. The thin dotted curve represents Pd calculated with the same potential 
but with modified r D (=1.55 fm). The thick (thin) dashed curve denotes Pp calculated with 
a F = 0.25 fm (a F = 0.45 fm). 

Fig. 3. The results of calculations for elastic probabilities, as a function of the reduced 
distance D for the 16 O+ 208 Pb system at E lab = 80, 84, 90, 96, and 102 MeV, using the final 
fully dispersive optical potential. 

Fig. 4. The Stelson plot of Si = ^E cm tji for direct reaction (i = D, solid circles) and fu- 
sion {i = F, open circles) cross sections. The straight lines are drawn to show the extraction 
of the threshold energies. Thin lines connecting the circles are only to guide the eyes. 
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Fig. 5. The strength parameters Vd(E) and Wd(E) for the direct reaction potential as 
functions of E. The open and solid circles are the values extracted from the x 2 -analyses. 
The solid lines denote Wd(E) and Vd(E) calculated, respectively, from Eq. (3.7) and from 
Eq. (3.5) together with Eq. (3.7). The thin lines connecting the circles are only to guide the 
eyes. 

Fig. 6. The strength parameters V F (E) and W F (E) for the fusion potential as functions 
of E. The open and solid circles are the values extracted from the x 2 -analyses. The solid 
lines denote Wf(E) and Vf(E) calculated, respectively, with Eq. (3.8) and Eq. (3.5). The 
thin lines are to guide the eyes. 

Fig. 7. The ratios of the elastic scattering cross sections to Rutherford cross sections 
calculated with our final optical potential for the 16 O+ 208 Pb system are shown in comparison 
with the experimental data. The data are taken from Refs. 3 and 4. 

Fig. 8. The direct reaction and fusion cross sections calculated with our final optical 
potential for the 16 O+ 208 Pb system are shown in comparison with the experimental data. 
The direct reaction data are taken from Refs. 3 and 4, while the fusion ones from Ref. 8. 
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